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In this paper we investigate certain properties of semi-totalistic cellular automata (CA) on the well 
known quasi-periodic kite and dart two dimensional tiling of the plane presented by Roger Penrose. 

We show that, despite the irregularity of the underlying grid, it is possible to devise a semi- 
totalistic CA capable of simulating any boolean circuit and any Turing machine on this aperiodic 
tiling. 



1 Introduction 

Since the discovery of quasi-periodic tilings and the related quasi-cristals, several researchers have stud- 
ied the behaviors of cellular automata (CA) on quasi-periodic tilings. 

In 2005, Chidyagwai and Reiter showed that the broken symmetry of quasi-periodic tilings, while 
still retaining a highly organized structure, could be used to simulate complex growth of snow crystals Q 
by studying real- valued CA on such lattices. They could produce global n-fo\A symmetry models where 
regular hexagonal grids could only produce 6-fold symmetry models. 

Reiter also investigated the behaviors of cyclic CA rules on several quasi-periodic tilings. Cyclic CA 
are special CA that can be regarded as a model of excitable media EE)). Once a cell enters a specific 
state, it must go through a fixed cycle of other states before being active again, thus implementing a 
refractory period mechanism. Reiter found that cyclic CA on quasi-periodic tilings, despite having a 
very simple local behavior, can take advantage of the lattice geometry to exhibit complex and smooth 
spiral patterns [14]. A variety of Cyclic CA rules are illustrated and commented in JT6l . 

Variants of John H. Conway's Game of Life [3] have also been studied on Penrose tilings. In particu- 
lar, Owens and Stepney found many stable and periodic patterns, similar to the ones known in the usual 
euclidian Game of Life |Q21[T3]|. Despite this, no equivalent to the glider structures of the euclidian case 
are known on quasi-periodic variants of the Game of Life. Indeed, because of the non-periodic nature of 
the underlying grid, the very definition of gliders is unclear. 

In this paper, we investigate some semi-totalistic CA on Penrose's kite and dart aperiodic tiling of 
the plane. We show in particular that there exists an 8-states semi-totalistic CA capable of simulating 
any boolean circuit. The usual difficulty of synchronizing signals, made even more problematic in an 
aperiodic setting, is avoided by embedding universal asynchronous logical elements. This construction 
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is similar to that employed to construct a universal CA on an hyperbolic space \7\. Once we are able 
to simulate logical circuits, we use Morita's construction JTOl to simulate a Turing machine with simple 
finite memory logical elements. 

2 Semi-totalistic Cellular Automata on Penrose Tilings 

2.1 Penrose's Kite and Dart Tiling 

This famous tiling of the plane was discovered by Roger Penrose as a simplification of his previous 6- 
tiles aperiodic tilings. This tiling is generated by two quadrilateral proto-tiles : a convex one called kite 
and a non-convex one called dart. To avoid periodic arrangements, tiling constraints are added on the 
tiles 0. 

This tiling has been widely studied, and its many properties will not be discussed in details in this 
article. 

There exist infinitely many different possible tilings of the plane with the kite and dart tiles. They are 
all aperiodic but quasi-periodic, that is to say that there exists a constant feEl such that for any r e R, 
any finite pattern of radius r that appears on the tiling appears in every area of radius kxr. 

2.2 The Moore Neighborhood on a Penrose Tiling 

In order to define a CA on a Penrose tiling, we need to define a neighboring relation between the different 
tiles. To do so, we extend the usual Moore neighborhood (8 closest neighbors on the regular square grid) 
by considering that two Penrose tiles are neighbors if they share at least one vertex. 

Because all tiles are not identical and not periodically arranged, different tiles might have different 
neighborhood. In the case of the kite and dart tiling, there are 8 different possible neighborhoods JT2l . 
illustrated by Figure [T] 

It is important to notice that not only these neighborhoods have different shapes, they also have 
different sizes ranging from a minimum of 9 cells to a maximum of 11. Although this would normally 
be a problem when defining a cellular automaton (because it would not be possible to define a unified 
rule for all cells), we will see that the special class of semi-totalistic CA is especially well suited for 
non-homogeneous neighborhoods. 




Figure 1 : The generalized Moore neighborhoods on a kite and dart Penrose tiling 111211 



K. Imai, T. Hatsuda, V. Poupet& K. Sato 



269 



2.3 Semi-totalistic Cellular Automata 

Because we want to work on aperiodic tilings of the plane, the usual definition of a cellular automaton 
must be changed in order to be compatible with the fact that different cells might have different neighbor- 
hoods. The subclass of semi-totalistic cellular automata is particularly well suited for this as the position 
of the different neighbors is irrelevant, only the number of neighbors in a given state has an effect on the 
evolution of a cell. 

Definition 1 (Semi-totalistic Cellular Automaton) A semi-totalistic cellular automaton is a pair $tf = 
(N, 8) where : 

• N EN is the number of states of the automaton. The set of states of the automaton is therefore 
£* = {0,1,...,N-1}; 

• 8 : x — )► is the totalistic local transition function of the automaton. 

Definition 2 (Configuration) A configuration of the automaton srf (N, 8) over a tiling of the plane by 
Penrose's tiles is a mapping £ that associates a state of to each tile. 

In accordance with usual cellular automaton terminology, in this context, the tiles of the Penrose 
tiling will be denoted as cells. 

Definition 3 (Global Evolution) Given a configuration (£, the image (£ f of £ by the automaton srf — 
(N, 8) is such that for each cell c, the state <£'(c) is defined by 

<£(c) = 8(q c ,no,ni,...,n N -i) 

where q c is the current state of c and for every i E [0, N — 1], n\ is the number of neighbors ofc in state i. 

We will denote by srf ((£) the image of a configuration £ by such a process. The global action of the 
automaton can be iterated over each new configuration, defining in such a way the infinite evolution of 
the initial configuration. 

In other words, at each time every cell of the automaton updates its state according to the number of 
neighbors it has in each state. 

Remark: In order to be as general as possible, and to be able to define a CA independently of the 
lattice on which we will consider it, we have defined the local transition function as a function from 
and infinite set £2^ x N^. However, once the geometric structure of the cells has been decided (in our 
case the Penrose tiling), it is sufficient to define the local transition function over the tuples that can 
actually appear. In particular, since we know that in the kite and dart tiling all cells have between 9 
and 11 neighbors (themselves included), only the input tuples of the form (#,no,ni, . . . where 
9 < < 1 1 should be considered. 

2.4 Previous Studies 

Owens and Stepney studied the Game of Life on Penrose tilings JT2l[T3l . 

Even though the average size of the neighborhoods is larger than that of the standard Moore neighbor- 
hood in the regular case, the same parameters for the rule (a cell survives with 2 or 3 live neighbors and 
a new cell is born when surrounded by 3 live neighbors) exhibits behaviors very similar to the classical 
Game of Life and many similar stable and periodic patterns were found. 

However, no glider or similarly propagating pattern has been found so far, possibly because of the 
strong irregularity of the encountered geometry as the signal moves. Gliders being the heart of the 
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computational and intrinsic universality proofs of the Game of Life in the classical case, this approach 
cannot be used in the quasi-periodic case. 

Another noteworthy difference between the Penrose case and the classical case is that long string-like 
(closed or not) connected alive cells can be stable on a Penrose tiling. Such patterns, called rings and 
snakes have been studied in |[T3ll . 



3 A Universal C A on Kite and Dart Penrose Tilings 

In this section we present the construction of a semi-totalistic CA working on Penrose's kite and dart 
tilings that is logically universal in the sense that it can simulate the behavior of any boolean circuit. 



3.1 Universal Reversible Logical Elements 

In order to construct a universal CA, we will use very simple reversible logical elements with 1-bit 
memory. 

The usual difficulty when embedding circuits in a larger structure (such as a cellular automaton) is 
to synchronize the different signals. In the case of an aperiodic geometry, the problem is even harder 
as construction patterns cannot be precisely shifted in a given direction. To solve this problem we use 
asynchronous elements that can simulate any circuit while only having one signal at a time. 

Figure [2] shows every 2- symbol reversible logical elements with 1-bit memory. There are only eight 
types of such elements. The left part of each illustration represents the element in state 0, the right part 
being state 1. Each element only accepts one signal from one of two possible inputs (left side) at a time. 
Lines (full and dotted) represent the output through which the signal will exit the element. 

When a signal traverses the element through a full line, the state of the element is changed. If the 
signal traverses through a dotted line its state is unchanged. 



2-0 



(eq. to wires) 



2^1 



(eq. to wires) 



2-2 




2^4 



2-5 (eq. to wires) 



2-16 (eq. to wires) 



2-17 



Figure 2: 2-symbol reversible logical elements with 1-bit memory ||9|. 



It has been shown that any reversible Turing machine can be constructed by 4- symbol reversible 
logical elements {rotary elements) [10|. Its computing process is executed by a signal, thus there is no 
need to synchronize signals. The rotary element can be realized by any combination of two elements 
amongst 2-3, 2-4, and 2-17 H1[TT]|. ^° ^ e combination of two elements are universal. 

In our construction we will use elements 2-3 and 2-17. 



3.2 General Overview of the Automaton 



The automaton has 8 states 
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• State is a quiescent state, meaning that a cell in state surrounded by other cells in the same state 
will not change. This state is therefore used as a "background" for the whole construction and a 
"filler" in the local transition rule to account for the different sizes of neighborhoods (the number 
of other states in the neighborhood is significant but the exact number of states is not). 

• States 1, 2 and 3 are used to transmit the signal through wire-like structures. Connected cells in 
state 1 will act as a wire along which the signal moves. The signal itself is made of two states to 
orient it. State 3 is the head of the signal and state 2 is the refractory state that remains after a cell 
was in state 3 in order to avoid the signal going back through the wire in the opposite direction. 

• Finally, states 4 to 7 are used internally in the logical elements to control the behavior and store 
the inner state of each element. 

Any boolean circuit can be simulated by an asynchronous 1 -signal circuit containing only copies of 
the elements 2-3 and 2-17, connected by wires HKTH. Because the circuit is two-dimensional, wires 
must also be able to cross. 

It is therefore enough to show that our automaton can accurately simulate the behavior of the logical 
elements and that these elements can be connected by wires. 

3.3 Simple Wires 

Wires can be easily represented by a continuous line of neighboring cells in state 1 surrounded by cells 
in state 0. It is important that each cell of the wire is connected to exactly two other cells of the same 
wire. This constraint is not very restrictive though, and it is easy to draw a wire going from any point 
of the tiling to any other. Simple wires cannot cross, and we will need a specific construction to allow 
crossings. 

A signal on the wire is represented by two consecutive cells in states 3 (front) and 2 (back). 
The transition rules needed to describe the behaviors on the wires are quite simple : 

• a cell in state 1 that has a neighbor in state 3 (head of the signal) changes to state 3 ; 

• a cell that is in state 3 changes to state 2 (refractory state) ; 

• a cell in state 2 returns to the simple wire state 1 ; 

• finally, cells in state 1 remain in that state in the absence of a neighboring signal. 
Figure [3] illustrates a wire on which there is a signal moving towards the right. 

3.4 Logical Elements 

The simulation of the elements 2-3 and 2-17 can be done by adding new states to the automaton. It 
would be very tedious to devise a construction that could behave as a specific logical gate in every pos- 
sible neighborhood, but we only need to have a pattern that works in a specific geometric neighborhood 
because the tilings produced by the kite and dart tiles are all quasi-periodic. 

This means that if we can position the widget at a given position because the surrounding tiles 
correspond to what we need, there are such appropriate positions in every sufficiently large surface of 
the tiling. 

For instance, Figures]?] and [5] illustrate a widget that behaves as element 2-3. This widget must be 
located on a neighborhood of type aO (according to the classification of Figure [T]) in order to work prop- 
erly as any other positioning would change the number of neighbors of one of the cells and significantly 
alter its behavior. 
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Figure 3: A signal on a wire and their state-color mapping table. 



But because of the quasi-periodicity of the tiling, there are such locations in every area of the tiling. 
We can therefore put different copies of the widget in different areas of the tiling and connect them with 
wires as needed (the wires do not require specific arrangements of the tiles). 

The quasi-periodicity can also be used for larger constructions. If we create part of a circuit in some 
specific location and want to copy it somewhere else, we know that there exists another appropriate 
location not too far from the original (the maximal distance is proportional to the radius of the widget we 
need to copy). Because of this, it is enough to verify that the widget works as intended in its appropriate 
environment to ensure that we can use it to build arbitrarily large circuits. 



state 



state 1 




a 
b 




Figure 4: Logical element 2-3 in states (left) and 1 (right). White cells in state 0, light gray cells in 
state 1. 



Figure [4] shows the two patterns that will simulate the element 2-3 with internal states and 1 respec- 
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tively. Figure [5] illustrates the evolution of this pattern when a signal enters from its second input (b) and 
exits through its first output (x) while changing its internal state. 

Similarly, Figures [6] and [7] illustrate the pattern used to simulate element 2-17. 





Figure 6: Logical element 2-17 in states (left) and 1 (right). White cells in state 0, light gray cells in 
state 1. 



Because we can control the environment in which the widget evolves, we can precisely set the local 
transitions up to behave as intended. 

Finally, in order to be able to connect all the elements freely, it is necessary to add the possibility 
of crossing wires. This is achieved by a special widget illustrated by Figure [8] and [9]. This last element 
is very similar to element 2-17 in its internal state 1, the only difference being that it should not change 
state when traversed by a signal. 

Table |3.4| gives a full description of the local transition rule of the automaton. 



3.5 Simulation of a Turing Machine 

It has been shown in lITOl that it is possible to create a circuit that simulates the behavior of a Turing 
machine cell by arranging eleven copies of a simple logical element with one-bit memory (called rotary 
element). In order to simulate a Turing machine it is then sufficient to connect an infinite line of such 
circuits. The resulting simulation only requires one active signal to move through the whole infinite 
circuit. 

The rotary element itself can be simulated by elements 2-3 and 2-17 but arranging an infinite line 
of copies of a given pattern cannot be done periodically on Penrose's tiling. However, since every tiling 
by Penrose tiles is quasi-periodic, if we can embed a circuit that simulates a Turing machine cell on a 
certain pattern of the tiling, there exists a constant N such that on every square portion of the plane of 
size N xN the pattern appears. This means that if we consider an infinite stripe of the plane divided in 
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Figure 7: Element 2-17: a transition from state 0, input a to state 1, output x. 




Figure 8: A crossing element. White cells in state 0, light gray cells in state 1. 
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Figure 9: A transtion of a crossing element. 
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Table 1: Local transition rule of the 8-states universal cellular automaton. 

The symbol is a wildcard that matches any value. In case of multiple matches, the first matching rule 
is applied (left to right, top to bottom). 
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squares of size N x N, we can put a Turing machine cell in each square and connect them into an infinite 



line as illustrated by Figure 10 




Figure 10: How to simulate the behavior of a Turing machine with copies of a circuit that simulates the 
behavior of one tape cell. 



4 Conclusion 

Despite working on an aperiodic structure with heterogeneous neighborhoods, the cellular automaton 
that we have described can simulate any logical circuit and even simulate the evolution of a Turing 
machine, if the infinite circuit is correctly embedded on the tiling (it can be, but in an aperiodic manner). 
It can work on any valid tiling of the plane by the kite and dart tiles because all such tilings contain the 
same finite patterns so any finite configuration on one can be inserted in another. 

The number of states of the automaton is likely not minimal for a universal CA. It seems likely that 
there exist ways to simulate the simple logical elements (or even a different universal combination) with 
fewer states but possibly larger widgets. It is however very difficult to optimize these constructions as 
many different environments have to be considered in the search. 

Intrinsic universality (the ability to simulate any other cellular automaton) should also be investi- 
gated. In that case however, a single-signal asynchronous system cannot work but it might be possible 
by arranging infinitely many circuits in similar environments, each simulating the behavior of a cell, and 
connecting them with wires of fixed length to avoid synchronization problems. 
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